Fast Computation of Bispectrum Features with Generalized Slow Roll 
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We develop a fast technique based on the generahzed slow-roll (GSR) approach for computing the 
curvature bispectrum of inflationary models with features. We show that all triangle configurations 
can be expressed in terms of three simple integrals over the inflationary background with typical 
accuracy of better than ~ 20%. With a flrst order GSR approach the typical accuracy can be 
improved to better than the ~ 5% level. We illustrate this technique with the step potential model 
that has been invoked to explain the WMAP temperature power spectrum glitches at f ~ 20 — 40 
and show that the maximum likelihood model falls short of observability by more than a factor of 
100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the 
local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum 
arises from integrals of nearly the same function of the background slow-roll parameters as the power 
spectrum but with a stronger weight to the epoch before horizon crossing. Ifence this technique 
enables reverse engineering of models with large bispectrum but small power spectrum features. 



I. INTRODUCTION 

In this paper, we develop the generahzed slow-roll 
(GSR) approach to obtain the bispectrum of curvature 
fluctuations produced by features in the inflaton poten- 
tial. Bispectra for these kinds of models have been pre- 
viously considered by a computationally intensive direct 
integration of the curvature fluctuations for each configu- 
ration [TJ[5]. The GSR approach provides a computation- 
ally efficient method that involves only a single function 
of the inflationary background. 

The GSR approximation was originally introduced to 
allow for accurate solutions of the power spectrum for 
models where the slow-rofl parameters are small but not 
necessarily constant [3]. This method was subsequently 
extended for cases in which the potential can have large 
features and the slow-roll parameters are not small |4j. 

Violations of slow roll arise when the inflaton traverses 
a feature, such as a step or a bump, in its potential. As 
the inflaton rolls across the feature, some of its potential 
energy is converted into kinetic energy or vice versa. By 
keeping the amplitude of the feature small, one is able 
to ensure that inflation is not interrupted; however, by 
arranging for the feature to be sharp, the inflaton un- 
dergoes a sharp transient acceleration which temporarily 
violates slow roll. 

The possibility of sharp features in the inflationary po- 
tential has a long history. Starobinsky first discussed the 
spectrum of adiabatic fluctuations for a potential with 
sharp features [5,. In particular, sudden downward step 
features can arise naturally in models of inflation de- 
rived from supergravity ?5 . Detailed numerical analysis 
showed that these features in the potential lead to oscil- 
lating features in the spectrum of curvature fluctuations 
[7] and consequently, such models have been invoked to 



explain glitches observed in the temperature power spec- 
trum at scales around £ 20 - 40 pfTT]. Although the 
bispectrum for such models has been previously studied 
[Tl[2], the intensive computation required has prevented 
a full assessment of its observability. The authors in [T] [5] 
wrote down an approximate analytic form for the bispec- 
trum produced by a step but did not pursue it further. 
Consequently we use this step model to illustrate the 
GSR bispectrum technique. While we find that the bis- 
pectrum produced by a step with the parameters chosen 
to best fit the glitch at £ = 20 — 40 is unobservable, by 
abandoning this prior an adjusting the width and height 
of the step, larger bispectra can be produced [T]. 

Features have also been shown to arise from other phe- 
nomenological processes: duality cascades [HI [13], water- 
fall transitions [14,, the imprints of heavy physics on the 
inflaton [T5] , fast phase transitions [TH] , and multiple field 
scenarios [17]. Sudden changes in the sound velocity in 
more general inflationary models have also been shown to 
give rise to such features [H] . The techniques we develop 
here apply only to cases where inflation is being driven 
by a single effective degree of freedom. The crucial as- 
sumption is that there is only one "clock." We work with 
canonical kinetic terms but expect the generalization to 
other forms to be straightforward. 

The outline of this paper is as follows. In ^jll] we re- 
view the formalism for computing the bispectrum in the 
literature. We derive the GSR approximation for the bis- 
pectrum in Appendi x [A] and apply it to the step model 
in §III| In Appendix |B[ we test our calculations against 
results in the literature. We use the fast GSR approx- 



imation in \ IV to estimate the observability of the bis- 
pectrum step features for a model that fits the WMAP 
power spectrum glitches. Throughout we work in units 
where the reduced Planck mass Afpi = (87rG)-i/2 = 1. 



II. BISPECTRUM FORMALISM 

In this section we begin by briefly outlining the general 
method used to evaluate the bispectrum employed by 
Maldacena [TS] and extended by Weinberg |^ . 

We work in comoving gauge, where the time slicing is 
chosen so that the scalar fluctuations are in the metric, 
and make use of the interaction picture, where these cur- 
vature fluctuations TZ evolve according to the equations 
of motion derived from the quadratic action 



dt(rxa^2eH 



(1) 



In this picture, the three-point function arises from inter- 
action terms defined to leading order by the cubic action 
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which evolve the states. Here 



<t)H 



2^2' 
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overdots are derivatives with respect to cosmic time t, di 
refers to the spatial derivative, spatial indices here are 
contracted with the Kronecker delta, — , and 

'. . .' denotes terms which vanish outside the horizon. 

The form of the action in Eq. ([2| is slightly differ- 
ent from that originally written down by Maldacena [12] . 
Maldacena omits the total derivative terms, while includ- 
ing terms proportional to the first order equations of mo- 
tion which are of the form f{Tl)SC2/STZ. However, since 
these terms are evaluated on-shell, their contribution to 
any Feynman graph is formally zero. Maldacena removes 
these terms by performing a field redefinition to a new 
variable TZn = TZ — f{TZ)- As described in [5T], the in- 
clusion of the boundary term in the last line of Eq. ([2| 
accounts for the terms one obtains by performing the field 
redefinition. 

We quantize the theory in the usual way. The transla- 
tional invariance of the background makes it convenient 
to expand the field TZ in Fourier components. 



7e(x,t) = 



d^q 

W) 



(4) 



Rotational invariance and Hermiticity imply that the 
most general solution takes the form 



7^q(^) = 7^,(^)a(q) + 7^:(^)at(-q) 



(5) 



where 9 = |q| is the magnitude of the comoving momen- 
tum and hats denote operators. In the quantum theory, 
a(q) and a^(q) can be thought of as creation and anni- 
hilation operators satisfying 



[a(q),at(q')] = (2^)35(q_q'). 



(6) 



With these definitions, we can construct the Fock space 
by applying creation operators to the state annihilated 
by all of the a(q) which we call the vacuum, |0). We can 
then define the two-point function, or propagator, 

(7^q(^)7^q,(^')> =7^,(^)7^;,(^')(2^)^<5(q-^q'), (7) 

which follows by simply normal ordering. The prescrip- 
tion of the "in-in" formalism for the expectation value 
of a product of field operators 0{t) is to evaluate the 
expression 



(0(t)) = {U\t,to)Oit)U{t,to)), 
where U{t,to) is the time evolution operator 



U{t,to) = Texp 
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The average in Eq. Q denoted by (...) is a quantum av- 
erage with respect to the vacuum state of the free field 
theory governed by the action in Eq. ([T]) . For the problem 
at hand, we take the initial time to to be in the asymp- 
totic past, to = — oo(l -|- ie), where the ie prescription 
projects out the Bunch-Davies state initially. 

The tree-level bispectrum is then given by expanding 
Eq. ^ [with 0{t) = 7^kl(^*)7^k2(^*)7^k3(^*)] to linear 
order 



(7^kl(^*)7^k.(^*)7^k3(^*)) 
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(10) 



In this work we are interested in potentials in which 
the inflaton undergoes a sharp transient acceleration but 
inflation is not interrupted. Consequently, e^f <C 1 ev- 
erywhere and the bispectrum is dominated by the term 
proportional to the derivative of rjH- To a very good 
approximation, for the purposes of this paper, the cubic 
action defined in Eq. ^ reduces to 



dtd X 



dt 



ahHieH-VH)n^n 



(11) 



For this cubic action in Eq. ([2|, the interaction Hamilto- 



nian is 



Hj{t) = 



d^3 



-j^(a'^eHieH-VH)n^n 
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Switching to conformal time, rj = Jt """^ dt' /a (defined to write this expression as 
be a positive quantity during inflation) and working in , 
Fourier space, we obtain for the interaction Hamiltonian 



Hiiv) 
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where here and throughout ' = d/dlnrj. In this expres- 
sion, the fields TZ are interaction picture fields whose time 
dependence is governed by the Hamiltonian derived from 
the quadratic action in Eq. ([ij. 

We are only interested in the connected part of the 
three-point function here, since this is the only part that 
contributes to the non-Gaussianity. This is evaluated 
from Eq. ( 10 1 using Eq. ([t]) and Wick's theorem. Defining 
the bispectrum through 

(7^k,7^k,7^k3) = (2TrfS{-k^+k2 + -k3)Bn{ki,k2,h), 

(14) 

we find 



(15) 



The bispectrum in Eq. ( 15 ) may appear to depend on 



the final time 77*. However, if we consider modes that are 
well outside the horizon, integrate Eq. ( 15 ) by parts and 



use the equation of motion for TZ, we obtain 



Bn{ki,k2, fcs) = m i{k{ + k^ + ki)Tlk,Tlk^'Rk, (16) 



We have also dropped the contributions from terms like 



dlwIZk 
dhir] 



(17) 



Bn{ki,k2,k^) =m{i{kl + kl 



kl) 



dT]a^eHieH - Vh) 
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which is obviously identically zero, regardless of the time 
at which we choose its evaluation. The badly divergent 
integral cannot acquire an imaginary part. 

One might worry that this statement holds only at 
leading order. In writing down Eq. ( 18 ) we have ne- 



glected the decaying mode which may be amplified by the 
divergent integrand if it were not decaying fast enough 
at late times. It is easy to see that this does not oc- 
cur by examining the asymptotic expansion of the mode 
functions TZ about the superhorizon limit [5D], 



7^,(r,) =TZl + 
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where TZ^ and Ak are complex constants and '. . .' refers 
to terms higher order in k in the expansion. Thus, the 
leading order corrections to Eq. ( 18 1 are proportional to 
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i[TZlAk] / di]a'^€H{eH - Vh) / 
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which, as long as the slow-roll parameters are well be- 
haved, converges essentially as an inverse power of the 
scale factor. 



Upon examination of Eq. ( 16 ), it can be seen that this 



conclusion is in fact a direct consequence of the gen- 
eral theorem proved by Weinberg [5D]. Eq. (16) can 



be thought of as arising from a "dangerous interaction." 
That is, an interaction that diverges as a at late times, 
but contains only fields and not time derivatives of fields. 
Then, due to Weinberg's theorem the integrals over the 
time coordinates of the interaction converge exponen- 
tially fast in cosmic time at late times. We also point 
out that the second term in Eq. ^ takes this form, and 
we are safe in concluding that it does not give large or 
divergent contributions to correlation functions at late 
times. 

This property shows that the time independence of the 
three-point function after horizon crossing is enforced by 
the boundary term regardless of how the slow-roll param- 
eters are evolving. 



since these converge essentially as inverse powers of the 
scale factor for these modes. Since the modes under con- 
sideration are well outside the horizon, they have become 
constant, and we can pull them out of the integral and 



III. GENERALIZED SLOW ROLL 

The generalized slow-roll (GSR) approach introduced 
by Stewart [3] is a technique for predicting the curvature 
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power spectrum in models where the inflation continues 
uninterrupted but the slow-roll parameters ch and 7]h 
evolve rapidly due to the presence of sharp features in 
the inflaton potential. 

Here we extend the GSR approach in two ways. First, 
we generalize the techniques for the calculation of the 
bispectrum. Second, we make them appropriate for the 
calculation of large amplitude features as was done for the 
power spectrum in Ref. [J]. The latter involves adding 
certain formally higher order terms, which can be justi- 
fied from an iterative expansion as discussed in Appendix 
\K\ We use this approach to develop a fast technique that 
can be used to approximate all configurations of the bis- 
pectrum for these kinds of models. 

We present only the main results in this section, leav- 
ing the details of the derivation and an overview of the 
GSR technique to Appendix |^ In the subsections we 
present comparisons of our approximation to the exact 
computation in the equilateral, squeezed and flat limits. 



A. Zeroth-Order Expressions 

The GSR approach proceeds by iteratively correcting 
the evolution of the mode function for the effect of devi- 



ations from de Sitter space [see Eq. (A7|]. 

At zeroth order, we employ only the de Sitter forms 
for the T^fc mode functions in a specific way described in 
Appendix [a] Eq. (15) then simplifies to 



BTz{ki,k2,k3) « -3-3-3 

1 2 3 

-Io{K)kik2k3-Ii{K)Y,kfkj 

+l2{K)K{kl + kl + ki)j (21) 
involving integrals separable in k, 

Jo V 

h{K) - GB(lnry,) + / ''^Gij(lnr,)cos(i^77), (22) 
/2(i^)^G.(ln,.)+r^^Gi,(ln,)^, 

which depend only on the perimeter of the triangle K — 
ki + k2 + k^. Thus the bispectra for all possible triangles 
can be efficiently obtained by precomputing these three 
integrals. Note that the trigonometric functions for Ii 
and I2 for Krj <^ 1 approach unity and so the expressions 
become independent of the arbitrarily chosen end point 
T]^. While this cancellation is guaranteed by the form of 
the action in Eq. ( 15 ), it is not guaranteed to occur order 



by order in the GSR expansion of the mode functions. 
This is due to the fact that in Eq. (15), the effects of the 



source are compensated by the response of the derivative 
of the curvature on super horizon scales. At zeroth order 



in the GSR approximation, the mode function does not 
respond to the feature at all, leaving this aspect of the 
source uncompensated. As we explain in Appendix [A] we 
have included the appropriate higher order terms in the 
source function in order to enforce this cancellation. 

The inclusion of higher order terms modifies the source 
from Eq. ([l5|) to 



where 



with G B 



V / 



/ = 



H 
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(23) 



(24) 



and 



k^P-ji/2TT^ is the curvature power spectrum. 
Note that for constant en ^ 1, aHrj = 1 and — A|j 
is the usual slow roll result for the power spectrum. 

The modification to the source and the replacement of 
the power spectrum for the zeroth-order external modes 



evaluated at 77* in Eq. ( 15 ) represent the two higher order 
corrections we have introduced. Both have the effect of 
enforcing that the bispectrum is insensitive to features 
in the inflaton potential for modes that are superhori- 
zon scale when the inflaton crosses the feature. For the 
source, this is achieved by making it a exact derivative of 
a combination of slow-roll parameters and affects mainly 
the Ii and I2 terms. 

In other words, we enforce that the three-point correla- 
tion function of curvature perturbations remains exactly 
constant outside the horizon at zeroth order in the GSR 
expansion of the mode functions. In Appendix |Aj we 
show that the higher order terms we add to the zeroth- 
order approach are compatible with a fully first order 
GSR computation. 

It is interesting to note the similarities and differences 
between how evolution in the slow-roll parameters affect 
the power spectrum and bispectrum. The source to the 
power spectrum is [see Eq. (A15)] 



Af f p 



(25) 



where the approximation follows for cases where /"// is 
the dominant contribution to both, as in the step poten- 
tial below. As noted above, in the ordinary slow-roll ap- 
proximation /^^ is the amplitude of the power spectrum. 
Thus, the presence of the extra factor of / serves to make 
the bispectrum scale as the square of the power spectrum, 
leaving an otherwise similar source for the power spec- 
trum and bispectrum. 

The main difference between the impact of features on 
the power spectrum and bispectrum is in the /q integral 
in Eq. (22 1. Unlike those for the power spectrum, this 
integral carries a divergent Kr]sm(Kr]) term as Kr] — > 00 
which can be traced back to the appearance of derivatives 



of the mode functions in Eq. ( 15 ) 
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FIG. 1. Bispcctrum source G'g and power spectrum source G' in 
the GSR approximation for the step potential model (see j jlll B\ . 
The bispectrum source has been rescaled by a constant related to 
Ag = 2.39 X 10~^ to eliminate its scaling with the power spectrum 
squared. Both sources share similar structure with features that 
integrate to zero as required for constant superhorizon behavior. 



FIG. 2. Bispectrum integrals in the zeroth-order GSR approxima- 
tion for the step potential model as a function of the perimeter of 
the triangles K = ki + k2 + ks. All bispectrum triangles can be 
formed efficiently from these three integrals. Weights reflect how 
the integrals contribute to equilateral configurations. Note that /q 
dominates the high k structure. 



The implication is that deviations generate bispectrum 
correlations while the modes are deeper within the hori- 
zon compared with the power spectrum. Thus one gener- 
ically expects that the impact of features in the inflaton 
potential on the bispectrum extends to higher k than in 
the power spectrum, enhancing their observability. 



B. Step Potential 



where As is an arbitrary constant that is of order the cur- 
vature power spectrum normalization without the step 
feature. In practice we take As = 2.39 x 10~^. Since G 
has dimensions of jWe typically divide it by some repre- 
sentative k^. In Fig. [2] we show the one-dimensional GSR 
bispectrum integrals of Eq. ( 22 ) for this model. Note that 



Io{K) dominates, especially at high k as discussed above. 



As an example, we consider the potential [7] 
1 



Vi<P) = -m 



2 j,2 



1 + c tanh 



(26) 



which corresponds to a smooth step ai (j) = (j)s of frac- 
tional height c and width d. Such a feature in the 
inflationary potential has been invoked to explain the 
'glitches' in the CMB temperature anisotropy data at 
£ = 20 - 40 pi [To]. Addressing the observability of the 
corresponding features in the bispectrum requires a fast 
approach to their calculation (see ^IV). 

For concreteness we adopt the WMAP5 maximum 
likelihood values of the parameters of the step po- 
tential dl] {m,c,d,(j),} = {7.126 x 10-^,1.505 x 
10-3,0.02705,14.668}. We plot the bispectrum and 
power spectrum source functions for this model in Fig. [T] 
Note their similar structure once rescaled in amplitude. 
In Appendix|B]we consider an alternate choice of param- 
eters to make contact with results in the literature. 

Following the notation of the existing literature [H [2] , 
we construct plots related to 



Q{ki,k2, ks 



"'l"'2'>'3 



(2^)4A| 



J-Bniki,k2,k3), (27) 



C. Equilateral Limit 



In the equilateral limit fci = k2 



kca and 



Eqs. (211 and (27) simplify to 



^(fcoq, fcoq, fccq) A^(fcoq)r fir i OTl 



1,3 



4A| 



(28) 



For equilateral triangles, the contribution from Iq dom- 
inates the result (see Fig. [2]) and hence the comparison 
with the numerical evaluation of Eq. ([Tsl ) tests this aspect 
of the GSR approximation (see Fig. Isl). The zeroth-order 
GSR approximation captures the amplitude and phase of 
the oscillations fairly well. The largest deviations are at 
low k for modes that were right on the horizon when the 
inflaton crossed the feature in Fig. [T] We show in Ap- 
pendix |A] that these are associated with the neglect of 
first order terms involving other combinations of the real 
and imaginary parts of the mode functions. 
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FIG. 3. Zcroth-order GSR approximation for equilateral configu- 
rations in the step potential model compared with exact results of 
evaluating Eq. 



D. Squeezed Limit 



FIG. 4. Zeroth-order GSR approximation for squeezed bispectra 
ki = kg <Si k2 ~ = kij in the st ep potential model compared 
with exact results of evaluating Eq. | |15| . Here fcg = 10~^ Mpc"'^. 
For comparison the prediction from the consistency relation with 
the slope of the curvature power spectrum is also shown. 



For squeezed triangles ks — ki <ti k2 ~ = and 



kl 



4i| 



[-2/1+4/2]^^ 



-.2k, 



(29) 



Note that in this hmit only Ii and I2 contribute and 
hence the comparison in Fig. |4] tests a different aspect of 
the GSR approximation. In particular Ii and I2 carry 
the main impact of the source correction discussed below 
Eq. (22) since their windows carry superhorizon weight. 
In Appendix |Xj we use these triangle configurations to 
develop and test our approximation. 

As a further check, it is well known that the bispec- 
trum of curvature fluctuations produced by an inflation- 
ary model with a single 'clock' obeys a consistency rela- 
tion which relates the squeezed limit of the bispectrum to 
the slope of the power spectrum [19] . The squeezed limit 
corresponds to one of the curvature fluctuations in the 
bispectrum having a much longer wavelength than the 
remaining two /eg <C fc^. In this limit, the consistency 
relation implies |22j 



kl 



din fc 



(30) 



In addition to the comparison of our approximation to 
the numerical results in this limit, we use the numeri- 
cally computed local slope to test the consistency rela- 
tion itself. The result is plotted in Fig. |4] which shows 
excellent agreement between the numerically computed 
bispectrum in the squeezed limit, and the result obtained 
from the slope of the power spectrum and the consistency 
relation. There are small discernible differences away 
from the effects of the feature; however, we attribute 
these to our truncation of the action at Eq. (11). 
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FIG. 5. Zeroth-order GSR approximation for flat configurations 
with kl = kp with k2 = k^ = kp/2 in the step potential model 
compared with exact results of evaluating Eq. (|15|. 



E. Flat Limit 

The final limit we consider is the flat limit, where kp = 
kl = 2k2 = 2fc3 and the wavevectors are co-linear. Here 



g{kF,kp/2,kF/2) _ AT^{kF)Al{kF/2) 



1,3 
ft, p 



16i| 



(31) 



X [-/o - 7/1 + l2l2]K^2k, 



and the approximation involves a combination of all three 
integrals. In Fig. [5j we show that the approximation 
works comparably well for flat triangles as equilateral 
and squeezed triangles. 
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IV. CMB SIGNAL-TO-NOISE 

Given the step potential model that fits the glitches 
in the CMB power spectrum, the corresponding features 
in the bispectrum described in the previous section 
[TT| are a firm prediction [U [5] . What is less clear is to 
what extent they are observable. Indeed that the model 
is designed to fit low multipole t ^ 20 — 40 glitches in 
the power spectrum implies that these features will be 
strongly impacted by the cosmic variance of the dominant 
Gaussian fluctuations [1]. On the other hand, bispectrum 
features generically extend to higher fc and hence I than 
power spectra (see pll A[ ). 

An estimate of the signal-to-noise in the bispectrum 
for the step potential has been hampered by the lack of a 
computationally efficient method for estimating the cur- 
vature bispectrum. Our zeroth-order GSR technique is 
ideal for these purposes as the bispectrum for any con- 
figuration can be simply formed from three precomputed 
integrals in Eq. ( 22 ) . 



A. Cosmic Variance 

The temperature or angular bispectrum is defined as 
the three-point function of the spherical harmonic coeffi- 
cients aira of the temperature anisotropy 



E 



^1 h ^3 

TOi m2 m3 



\rn\ 



(32) 

The cosmic variance of the Gaussian part of the field puts 
an irreducible limit on the signal-to-noise ratio of 



E 



(33) 



where 

^111213 = [1 + ^ili2 + ^(2^3 + ^t^ll + "^^1112^1213] (34) 

accounts for permuted contractions of repeated ^'s and 
the angular power spectrum is defined by 



(35) 



Thus to evaluate the signal-to-noise in the angular bis- 
pectrum we require not only a fast method for computing 
the curvature bispectrum but also for computing angular 
bispectra from curvature bispectra. 



In the flat-sky approximation, the angular bispectrum 
is defined by the three-point function of the Fourier mo- 
ments of the temperature field given by a(l) 

{a{h)a{h)a{h)) = (2^)2^(li +h + l3)%„,„,3). (36) 

For £1, £2, ^3 ^ I7 it is related to the all-sky bispectra as 



Bf 



(2£i + l)(24 + l)(2£3 + l) ( £2 4 





An 



xB, 



(37) 



Under the Sachs-Wolfe approximation, the tempera- 
ture field as a function of angle n on the sky is the pro- 
jection of the curvature field onto the sphere at the re- 
combination distance D from the observer 



a(n) ^ --7^(x = Dn) = 

5 5 



ik-Df 



, (38) 



so that 



dh 



'n 



where 



53L»4 27r }_ 

Y.Bn{k\MM) 



ki - (li/i5,fci||), 

k2 - (l2/i^,fc2||), 

ks = -ki - k2, 



27r 



(39) 



(40) 



and II is the direction along the line of sight, orthogonal 
to the plane of the sky. Note that even in this approxi- 
mation, the signal-to-noise in the angular bispectrum is 
a five dimensional sum over the curvature bispectrum. 

For consistency, we also compute Cg under the same 
fiat-sky, Sachs- Wolfe approximation for the cosmic vari- 



ance in Eq. ( 33 ) 



52^)2 



2^ 



P7^(k=(l/Afc||)). 



(41) 



In Fig. § we show the cumulative signal-to-noise as a 
function of the maximum I in the sum of Eq. ( 33 ) for the 
step model of §111 B| This model falls short of predicting 
observable effects by more than 10^ in the number of 
bispectrum triangles or more than 10^ in the amplitude 
of the bispectrum. Consequently, the crudeness of our 
curvature to angular bispectrum calculation is justified. 



B. Approximations 

To obtain an order of magnitude estimate for the 
signal-to-noise ratio, we seek only a crude computation of 
the angular bispectrum from the curvature bispectrum. 
We therefore take the flat-sky approach and the Sachs- 
Wolfe limit for the temperature anisotropy. 



V. DISCUSSION 

We have developed a computationally efficient tech- 
nique based on the generalized slow-roll approach for 
calculating the curvature bispectrum of models with fea- 
tures in the inflaton potential. This technique allows 
all configurations of the bispectrum to be calculated 
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FIG. 6. Approximate cumulative signal-to-noisc in the CMB an- 
gular bispectrum for the step model assuming cosmic variance lim- 
ited temperature measurements out to a maximum multipole ^max- 
This model falls short of observable by more than 10^ in the number 
of bispectrum triangles or 10^ in the amplitude of the bispectrum. 



based on three precomputed integrals over the inflation- 
ary background reducing the dimensionaUty of the prob- 
lem from three to one. 

In this zeroth-order approximation, the accuracy is suf- 
ficient to capture the overall amplitude and structure of 
the bispectra to typically better than 20%. We have also 
developed a first order approximation that that brings 
the typical accuracy to better than several percent at 
the expense of raising the dimensionality to two. 

The accuracy of the zeroth-order approximation more 
than suffices to make an estimate of the observability 
of the bispectrum features in the step potential model 
that best fits the WMAP power spectrum glitches at 
£ ^ 20 — 40. We find that the bispectrum amplitude 
is more that a factor of 100 too small to be observable in 
a cosmic variance limited measurement of CMB temper- 
ature anisotropy. 

We have also explicitly verified that the consistency 
relation holds for a model of inflation which violates slow 
roll due to a sharp downward step feature in its poten- 
tial. We find excellent agreement through the region af- 
fected by the feature, while small disparities at smaller 
and larger wavenumber are consistent with our neglect 
of additional terms suppressed by 0(e/f ) in the cubic ac- 
tion. 

Two key general insights arise from the comparison of 
the GSR approach to the power spectrum and bispec- 
trum. First, the dominant source of both is a similar 
combination of slow-roll parameters from the solution of 
the inflationary background. Second, the bispectrum is 
generated while modes were deeper within the horizon 
compared with the power spectrum. 

The latter fact implies that features in the bispec- 
trum generically persist to higher wavenumber k than the 
power spectrum and provides a template for constructing 
models with large bispectrum features but small power 



spectrum features. 

Indeed recent work has revealed an intriguing 'decou- 
pling' limit under which large non-Gaussianity may be 
produced while the power spectrum remains largely un- 
perturbed [m [5S] . The application of the techniques we 
develop here may be able to extend the region of validity 
for these scenarios into the region where gravity is not 
completely decoupled, merely very weakly coupled. We 
defer these considerations to a future work. 
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Appendix A: Generalized Slow Roll 

We begin in |A l| by reviewing the generalized slow-roll 
(GSR) formalism collecting some important re- 

sults and establishing our notation. Employing the same 
methods used to correct the power spectrum in lA 2| we 



derive the GSR expansion in for the bispectrum in |A^ 
We apply this expansion to the zeroth-order expansion 
used for the main results of the paper in |A4| This ex- 
pression involves corrections for the superhorizon modes 
based on a more accurate first order approach that we 



develop and test in A 5 



The treatment and notations used here mirror the pre- 
sentation in [3]. The reader is directed to these earlier 
works for further details. For notational compactness we 
express functions of the three separate bispectrum k val- 
ues as 



F{ki) = F{k^,k2,k:i) 
throughout this Appendix. 



(Al) 



1. GSR Corrections to the Mode Functions 

Varying the action ([!]), the mode functions, 7lk{rj), sat- 
isfy the Mukhanov-Sasaki equation 



l-^)^y^+n,^0, (A2) 
J J X ax 



where x = krj and / were defined in Eq. ( 24 1 
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We seek to expand the mode functions in a perturba- 
tive series around their infinitely slow-roh or de Sitter 
form. Since the comoving curvature TZk is undefined in 
that limit, it is useful to change variables to 



27r2 X 



(A3) 



This variable yi also carries the interpretation of the in- 
flaton field fluctuation in the spatially flat gauge and is 
well defined in the de Sitter limit. 

With this correspondence it is easy to see that well 
inside the horizon yi is more immune to features in the 
inflaton potential than TZk as it represents a free field. 
The converse is true outside the horizon. The latter is the 
primary flaw in GSR that we seek to rectify by enforcing 
constant curvature outside the horizon. 

As an aside, note that similar superhorizon issues arise 
even for the computation of the bispectrum with the ex- 
act TZk- Although TZk does not respond significantly to 
the feature on such scales, TZ'^, does and in such a way to 
exactly cancel the behavior in the £h, t]h sources to keep 
the bispectrum constant in Eq. (f5). 

We then have 



where 



1 
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r - if 

f ' 



(A4) 



(A5) 



and primes denote derivatives with respect to In 77. 

The homogeneous Eq. (A4) corresponds to exact de 
Sitter space {H = const.) and has solutions 



(A6) 



yo{x) = 1 + 



and ?/q (x) which depend in the same way on a; = kj] for all 
k. Given these solutions, we can use the Green function 



of the homogeneous operator to invert Eq. ( A4 1 , 

yi{x) = yo{x) + L{x, u)yi{u), 

with 



(A7) 



L{x,u)yi{u) 



du 



g{\nf])y,{u)'^[yl{u)yo{x)], 



(A8) 



where u = kf}. Note that unlike j/Qj Vi is not the same 
function oi x — k-q for all k. 

Then, assuming that the new solution is not too dif- 
ferent from the de Sitter space result, we can employ the 
Born approximation and to solve the formal solution in 
Eq. IATI) iteratively 



yi{x) = yo{x) + Lix,u)yo{u) + L{x,w)L{w,u)yo{u) + ... 

(A9) 



We will also use the notation 

3, 



Wiu) 



-^[yoiu)]?Si[yoiu)] 
a 

3sin(2u) 3cos(27i) 3sin(27i) 



2w3 



2u 



X{u) ^-^[yo{umyo{u)] 



3cos(2w) 3sin(2u) 3cos(2w) 



and 



2^3 



1 u 

i[yo{u)Myf,{u)] = 1 + ^ - -X{u). 



2u 



(AlO) 



(All) 



In the limit of small u, these window functions behave 
as lim^^o W{u) = 1 and lim„_j.o A^(w) = u'^/3. 

2. Power Spectrum Expansion 

The curvature power spectrum is given by 

A'^{k)^\\m^y,{x)y*{x). (A12) 

At zeroth-order in the mode function correction and first 
order in g we obtain [3] 



lim (xyO =i - ^ 



+ 



u w 
du 



u 
du 
u 



^gilnrf) 
W{u)g{\nr]) 



(A13) 



A(u)g(ln7?)+0(x2), 



(which corrects a sign typo for the third term in Eq. (19) 
4J) where x — krj^, and u — krj and so 



1 



2f 
3 / 



drj 
V 



W{kri)g(\nr]) 



(AM) 

since the X term contributes quadratically. 

The problem with this expression is that the power 
spectrum depends on the arbitrary end point of the in- 
tegration r|.^,. The origin of this problem is that the cur- 
vature computed in the GSR approximation is not guar- 
anteed to be constant outside the horizon and constancy 
is only enforced order by order in g [4 . With large fluc- 
tuations, contributions that are formally higher order in 
g supply necessary corrections. 

To find these corrections, we can compute all of the 
first-order mode function corrections [26j and retain the 
ones that contribute on superhorizon scales [1] 



InA^SR(fc) = G(lnr;, 



drj 



W{kTj)G'{lmj), (A15) 
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where the modified source is 



2 / f'^ 



(A16) 



Formally G' involves a first-order correction to the source 
since /'// — C'(g). Note that for superhorizon modes 
when the inflaton crosses the feature W{k'ri) — > 1 and 
the power spectrum no longer depends on the arbitrary 
end point. 

We will follow this procedure to define our zeroth-order 



bispectrum formulation in Eq. (21 1 



3. Bispectrum Expansion 

In terms of the above notation, the bispectrum in 
Eq. ( 15 ) takes the form 



+ Boundary Terms 



-gsi 



J (In ry)D^ [yl {kiri)y*2 {k2ri)yl [k^^r])] | 

(A17) 

and the differential operator 
d 



Dr. 



din J] 



3 1 



11. 
f 



(A18) 



In the Eq. ( A17), "Boundary Terms" refers to the second 
term in Eq. (15). While these terms are physically im- 



portant, they represent a small correction and we defer 
their discussion to a subsection. 

Here the unmodified source function, which we will 
subsequently modify in a similar prescription to g ~> 
3G72, is 



(eg - flHY 
f 



(A19) 



By analogy to the power spectrum window functions, let 
us define the window functions 



Wb{xi) = di[yi{xi)y2{x2)y3{x3)], 
Xb{xi) = ^[yi{xi)y2{x2)y3{x3)]- 



(A20) 



For the zeroth-order mode functions, yi — >■ j/Oi W^s ~^ 
Xb -> Xl 



W%{x,) 



mx^) 



1 



{ — X -\- XiX2X^) COS X 



XlX2X^ 

+{1 — X2X3 — X1X2 — 2^12:3) sinX 
1 

X1X2X3 

+{X — X1X2X3) sinX 



(1 — X2X3 — X1X2 — X1X3) cosX 

(A21) 



where X = xi + X2 + x^, and is not to be confused with 
the window function X{kri). Note that the trigonometric 



functions depend only on the perimeter of the triangle K 
and that the prefactors are pure powers of fc^. This is 
a critical simplification achieved by using zeroth-order 
mode functions. 

As we shall see, their superhorizon limit will be par- 
ticularly important 



7,3 I 7,3 



k^^Q " 3kik2k3 

1 



(A22) 



'-+0{kfr^% 

fc.-i-o fciK2K3?7 2 kik2k^'q 

To an arbitrary order, the bispectrum can be written 

(2^)' Vl 



4fclfc2fc3 P 

drj 



^{i[WB{h7u) + iXB{hT^,)] 



V 



gB{lnv)D^ [WBihr,) - iXB{k,r^)] }. 

(A23) 

Specifically, using this notation we can write the zeroth 
and first-order bispectrum as 

(27r)4 ^3 .00 



dr] 



gB{\^rj)D,j{ 



(A24) 



Akik2k3 p J^^ 77 

- X%{k,r],)W°ihi^){l - 5R[C(A;„77) - C(fc„ ??*)]) 
+ WUkiV*)WUhv)^[C{h,^,)^C{h,r,)] 
+ X%{hrj,)X%{h7j) 3[C(fc„ ry*) - C{h, 77)]}, 



where 



C{h,r]) = Yl 



djigjlnfj) yoikjfj)^^ , 



^[yoik3V)yo{kjri)]- 
(A25) 

The fact that C depends on the three fc's presents the 
main obstacle to simplifying the first-order expressions. 



4. Zeroth-Order Bispectrum 

We define the order of our approximation by whether 
the mode functions used in the construction are the 
zeroth-order de Sitter yg or an iterative expansion. Thus 
we allow ourselves the freedom to add higher order cor- 
rections to the source function but not the mode function. 
This operational definition, rather than order counting 
in g is motivated by simple form of the bispectrum that 
results. With only zeroth-order mode functions, all tri- 
angles can be expressed in terms of single integrals that 
depend only on the perimeter of the triangle. 



a. Unmodified Source 

We begin with the zeroth-order approximation using 
the original source gB in Eq. (A19). This form is consis- 



tently first order in the deviations from slow roll and of 
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I I 

Unmodified GSRO 

Exact 




1 



0.001 

k, (Mpc-1) 



Squeezed _ 

0.01 



FIG. 7. Unmodified source with zeroth-order GSR approximation 
for squeezed triangles as in Fig. [4] SoHd fines denote unmodified 
form ( |A26[ |. Red dashed fines represent the replacement of the 
externaTl// with A-^^ in Eq. ||A28||. Here rjt = IMpc. 



the terms in Eq. (A24| only the X%{kiri^,)Wg{kiri) lead- 
ing order term survives 



dr] 



gB{\nrj){W'^' + 3W^>) 



(A26) 



In Fig. [7j we compare this approximation with the exact 
result for squeezed triangles. There are two obvious flaws 
in this expression. First the external / = /(In 77*) factors 
depend on the arbitrary end point. Second, in the super - 
horizon limit Wg goes to a constant defined in Eq. ( A22 1. 



Thus, analogous to the behavior of the power spectrum 
in the same approximation, the feature is imprinted on 
superhorizon modes. 



b. Source Modification 

Just as in the power spectrum case, we can fix these 
problems by examining the first-order mode function cor- 
rections in the superhorizon limit. We start from the full 
first-order expression Eq. ( A24 ) and extract a factor of 



l~di[C{k,,rj,)] = l + ^+Y, 



W{k,v), 



(A27) 

where we have evaluated the expression in the limit 
kiTj^, ^ 1. We have also dropped higher order terms. 

Combined with the factors of 1// we recognize these 
factors as the GSR expression ( A14) for the power spec- 



trum. Our prescription therefore is to replace these terms 
with the exact power spectrum 



This fixes the problem of the external / factors in 
Eq. ( A26 1 . In Fig. [jj we show the impact of applying 
this correction. The low k features in particular are mod- 
ulated and enhanced by the external power spectra but 
the superhorizon problem remains. 

We then look at the superhorizon limit of the rest of 
the first-order corrections. Since Wb scales as a constant 
in the superhorizon limit, we need only keep the constant 
parts of the expansion of '^[C{ki^ rf)] 



lim mC{h,ri)] 



3 

E 

1 



— -^g{\xi-q) 
rj rj-^ 



(A29) 



'^W{k,fi)g(\ni)) + Oik^.ff) 
rj ■' 



On the other hand, Xb behaves as k^ so we need to 
keep up to 0{kf) parts of 3[C(fci, 77)]: 



lim 5[C(fci, 77)1 



2{k,rjf 



{kjrif 



00 



r, V T 
5^2 --3\ 



'^X{kjf])g{\afi) 
V 

^ dfi 

4wik,fj)gilnfi) 
V 

dfjr]^ 

-—g(\nr]) 



+0{k'^T]'fl 

Then, working to order 0{gBf'/f) 
3(27r^4 



(A30) 



drj 



An{ki)An{k2)Anik3) 



kl 



-9b 



Now, to a good approximation. 



9B = 



f 



fl f 



1 [f 
f\f 



3fcifc2fc3 

(A31) 



(A32) 



With this approximation for the source, integrating the 
second line of Eq. (A31) by parts one obtains 



V \f f 



1/ r 

f f Jrj, V f 



(A33) 



Dropping the slow-roll suppressed second term, we can 
account for this first-order effect on superhorizon scales 
by the replacement 



9b 



1 [f 

f \ f 



If 
f f 



eg - Vh 
f 



— G'b- 



1 



(l-3?[C(fc„77,)]) = AK(fci)AK(fc2)AK(fc3). (A28) term in Eq. ([15 



(A34) 

In the superhorizon limit, the contribution of the first 
to the bispectrum then reduces to the in- 
tegral of a total derivative of slow-roll parameters, which 
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are small at late times, once the inflaton has settled back 
onto the slow-roll attractor. 

To summarize, our zeroth-order approximation con- 
sists of the replacement of the source by Eq. (A34) and 



the replacement of the external factors using Eq. (A28) 



while using the de Sitter mode functions everywhere else. 
These replacements lead to the expression in Eq. (21) 



with the help of the explicit form for the Wg window 
in Eq. (A21|. In particular note that the Iq term arises 



from the W[ 



term in Eq. ( A26 ) 



Boun 



Terms 



While small, the boundary term in Eq. ( 15 ) hitherto 



omitted plays an important role in ensuring that the bis- 
pectrum becomes independent of time on superhorizon 
scales. 

The contribution from the boundary term can be ob- 
tained from the results of the previous section with the 
replacement 



drj 
V 



5B(ln??) 



(eh - Vh) 



f 



(A35) 



together the replacement of all 77 by 77,^. Using this re- 
placement on Eq. (A31) yields for the boundary term 



7,3 



eg - Vh 
f 



1 - 



V f 



3fcifc2/C3 

(A36) 



Notice that, when added to Eq. (A31| using Eq. (A33l, 



one obtains perfect cancellation leaving a time indepen- 
dent result. Since the integral term in Eq. (A36) exactly 



cancels the second term in ( A33 ) we choose to omit both 
in practice. This brings the net boundary term to the 
G'B(ln?7*) contribution to Eq. (21). 



We also point out here that, while for convenience we 
have evaluated our expressions at a time when the infla- 
ton has passed the feature and is back on its slow-roll 
attractor, our results are not restricted to this regime. 
As we pointed out in Section |ll] the cancellation we are 
enforcing is exact and thus truly independent of the time 
at which the bispectrum is evaluated, even if this time is 
chosen to be when the slow-roll parameters are not small. 



5. First-Order Bispectrum 

We can further compute to first order in the mode 
function correction on all scales. Note that the first or- 
der corrections represent a calculation of the bispectrum 
to second order in the slow-roll parameters. Additionally, 
note that the boundary terms are always suppressed by 



slow-roll parameters evaluated at times well after the in- 
fiaton has traversed the feature which means that they 
are computationally irrelevant. 

Neglecting boundary terms we obtain at first order in 
the GSR approximation, the full bispectrum 



(2^)4 AGSR(fcl) AGSR(fc») AGSR(fc3 



^2 



drj 



+ UiD + UiE\{hri), (A37) 



where 
UiAihri) 



din?/ 

-j^ poo 



3 ) 3fi[2/o(fciry)yo(fc2'7)2/o(^3»y)], 
d 



"^^C (In fi)X{hfi) 



UisikiV) = - 



X '^[yo{kiv)yoik2'n)[yoih'n) + yoihv)]] 

-I- eye, 
2 J, J V 

-j^^ +3^ n[yo{kii])yo{k2v)yo{k3r])] 
+ cyc., 

Uic{hv)=-l r ^G'{\nfj)Xik,fj) 

X (^1^ + 3) ^yo{kiv)yo{k2v)yoik3v)] 

+ cyc., 
3 

4 

d 



UiDihri) 



1 



1 



^G'(lnfy) 1 , ~ i 
7, V yksV [ksV) 



UiEihri) 



dlnr] ^[2^o('^i^)2'o('^2'7) 

[Voikav) + ^0(^3??)]] + eye, 
- i^[yo{kii])yo{k2T])yo{k3V)] 



f 



1 - 



1 



■I 

^9B \f 



(A38) 



where eye. denotes the 2 additional cyclic permutations 
of the k indices. 

In order to fix a small second order correction to the 
superhorizon results we have replaced g — >■ 3G"/2 in the 
the first order expressions and likewise replace the source 
term 



f I 



±^ [f 
29B \ f 



(A39) 



in UiE- Furthermore for consistency, we use the GSR 
power spectrum approximation (A15) for the external 
terms. 
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1 
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11/ 
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FIG. 8. First order GSR mode function approximation for the 
step model compared with the exact solution. 



FIG. 9. Bispectrum in the step model of Ref. [1] [2] and Ag = 
4.415 X 10"^'' (see text). Our approximations work equally well for 
models with larger bispectrum features. 



For the case of the equilateral configurations, this ex- 
pression simplifies considerably and are practical to eval- 
uate 

(A40) 



where 



Uq^{v) — — 6cos3u + 



3fl-'y2)sin3u 



Ul\{v) = --[{?, + v'^)cosv + ?,{v'^ ~ l)cos3i; 



2v 



-|-2w(sin u — 3 sin 3f )] 

X / ?G'(lnry)X(fcr)), 
Jv. ^ 



(3 + v^) sinv 



2 cost! 



X / ^G"(ln77)VK(fcry), 

- — h V COS V + 2 sm V 
2 [\v J 

X / ?G'(lnry)X(A:77), 
9 

— |"(3 -I- w^) COS V + 3(^2 - 1) cos iv 

4:V 

-|-2w(sin w — 3 sin 3f )] 
3 

-3)cos3i; + (1 -3u2)sin3w] 



/' 



1 - 



(A41) 



We show the result of computing these terms for the 
step model of §111 B| in Fig. [8] Note that most of the er- 
ror in the zeroth-order approximation is corrected by the 
first order expression. The dominant correction is from 
U'^'X and involves the product of two independent inte- 
grals. Since its computation is no more intensive than 
the zeroth-order expression, it may be simultaneously 
computed as a monitor of the accuracy of the zeroth- 
order expression. Note that the integral of G'X is the 
same as the one that monitors the accuracy of the power 
spectrum approximation |?7j . 

This first order mode function correction breaks the 
property that all triangles depend only on integrals in- 
volving the perimeter K. Thus though more accurate it 
is of more limited utility for fast computations in that the 
six integrals involved are a function of two variables: the 
perimeter K and the wavenumber of the corrected mode 
function ^3. Nonetheless it still reduces the complexity 
of the exact computation which requires all three ki as 
well as the exact mode functions for each individual fc, . 



Appendix B: Comparison to Previous Work 

Previous calculations of the step bispectrum employed 
the parameters of: {m,c,d,(j)s} = {3 x 10^^,1.8 x 
10-3,0.022, 14.84} There are three importa nt dif- 

ferences between this model and the one defined in §111 B| 

First due to an error in the setting of the initial con- 
ditions in Ref. [H [10] the location of the feature 4>s is 
shifted and the corresponding features in the power spec- 
trum and bispectrum appear at lower k. Secondly, the 
maximum likelihood fit to WMAP5 data prefer a slightly 
larger width d compared with WMAP3 data. Since the 
bispectrum is particularly sensitive to the width d, its 
features appear more prominent and persist to higher k. 
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Finally Ref. [Hll] arbitrarily set the parameter m whereas 
Ref. [21 Ho] normalized to WMAP3. As a consequence, 
we take As = 4.415 x 10^^" to be compatible with this 



choice. 

In Fig. [9] we show our zeroth-order, first-order and 
exact calculation of this modeL We have verified that 
the exact calculation accurately reproduces the results 
in Ref. [T]. 
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